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Abstract 

We construct the effective potential for a QCD-like theory using the auxiliary field 
method. The chiral phase transition exhibited by the model at finite temperature and 
the quark chemical potential is studied from the viewpoint of the shape change of the 
potential near the critical point. We further generalize the effective potential so as 
to have quark number and scalar quark densities as independent variables near the 
tri-critical point. 
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§1. Introduction 



The phase structure of quantum chromodynamics (QCD) at finite temperature T and 
quark chemical potential \x has has been investigated since the invent of QCD. 1 ) Today 
collider experiments using ultra-relativistic heavy ion beams are creating highly-excited QCD 
matter in the laboratory. 2 ) Extremely dense, cold hadronic matter is relevant to the physics 
of the inner structure of neutron stars. 3 ) 

It is currently an accepted concept that the chiral symmetry of QCD, which is spon- 
taneously broken in the vacuum, will be restored at sufficiently high temperature and/or 
quark chemical potential. Although steady progress has been made in the study of lattice 
QCD at finite T and /i, numerical simulations are still hard to be carried out with physically 
realistic conditions, especially for finite /i. 4 ) Many model studies have been undertaken and 
have provided information concerning the in-medium properties of QCD. 5 ) In Ref. 6), the 
T-/i phase diagram was investigated with the Nambu-Jona-Lasinio model, and evidence was 
found suggesting the existence of the tri-critical point. The phase diagram was also stud- 
ied using the Schwinger-Dyson (SD) equation for the fermion propagator with a variational 
ansatz in a QCD-like theory. 7 )' 8 ) With the possibility of color superconductivity being newly 
suggested in Refs. 9)— 12), the phase structure at high density but low temperature has been 
re-investigated by using the diquark condensate as a new order parameter, in addition to 
the usual quark-anti-quark qq condensate. 

The effective potential is a useful tool, and it provides a clear picture regarding the 
nature of the phase transitions in these studies. In this paper we re-investigate the chiral 
phase transition of a QCD-like theory, focusing on the shape change of the effective potential 
near the critical point. 

The QCD-like theory 13 H 5 ) is the renormalization-group (RG) improved ladder approxi- 
mation for the SD equation of QCD. This model describes the dynamical chiral symmetry 
breaking in the vacuum while retaining the correct high energy behavior of the quark mass 
function. Using this model and the low-density expansion, the pion-nucleon sigma term and 
the quark condensate at finite density have been calculated. 16 ) The results are consistent 
with those of other models. The chiral phase transition at finite T and/or /x in this model was 
studied in Refs. 17)-20). The color superconducting phase has also been investigated. 21 )' 22 ) 

The Cornwall- Jackiw-Tomboulis (CJT) potential functional 23 ) has been used in several 
works 8 )' 19 ) attempting to determine the non-perturbative extremum solution for the model. 
For instance, in Ref. 19), the mass function was assumed to have a certain trial form char- 
acterized by the value of the scalar qq condensate, which was determined by extremizing the 
CJT potential variationally. However, the interpretation of the CJT potential away from the 
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extremum is not obvious. 24 ** The solution is known to be a saddle point of the CJT action 
with respect to general variations. 

Instead of following the approach outlined above, we construct an effective potential using 
the auxiliary field method in the QCD-like theory. Within the stationary phase approxima- 
tion, the solution of this potential is known to be a local minimum, 24 ) and it coincides with 
the CJT potential when the external field is turned off. This is a suitable property for a 
variational method, and it could be applied to other order parameters, such as the diquark 
condensate, provided that the gauge invariance is treated properly. 

This paper is organized as follows. In the next section, we review the QCD-like theory in 
the vacuum and at finite T and /i. In §3, we define the effective potential for the qq condensate 
in the QCD-like theory. Furthermore, we generalize the potential such that it possesses the 
quark number density as the second order parameter, in addition to the qq condensate. 
We investigate the transition points through consideration of the global properties of these 
potentials. The last section is devoted to a summary. In Appendix A, we demonstrate the 
validity of our definition of the effective potential in the Nambu-Jona-Lasinio model. 

§2. The QCD-like theory 

In this section, we briefly review the QCD-like theory, 13 H 5 ) introducing the notation and 
approximations used in this paper. 

2.1. In the vacuum 

Let us start by deriving the SD equation in the QCD-like theory as an extremum condition 
for the effective potential obtained using the auxiliary field method at zero temperature and 
zero chemical potential. 13 )~ 15 ) 

First, we integrate out the gluon field in the expression for the partition function Z, 
omitting the gluon self-interaction term in the QCD Lagrangian with massless quarks: 

Z oc J Vi^V^VAexpi J d 4 xC 

ip(p)p i tp(p) 



VipVip exp 



v 



\ J (p - |) % (p + |) K^' s (p, k)^ {k + §)*,(*-§ 



(2-1) 



where ip{p) is the quark field in the momentum space and f p = f d A p/{2n) A . The indices, 
a, (3, ■ ■ ■ correspond to the Dirac structure, and T a is the color su(N c ) generator. Here, we 
defined the kernel K as 

K a ^ 5 (p, k) = g 2 (^T a ) 5a ( lu T a )^iD^(p - k) , (2-2) 
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with the gluon propagator iD^ v {p) given by 

iD^ip) = 9 - { , >P P /P . (2-3) 

The kernel K represents the gluon exchange between the quarks. Hereafter, we employ the 
Landau gauge (i.e., a = 0). The non-Abelian nature of the gluon interaction is treated in this 
model as the one- loop running of g in (|2-2|) : g(max(p E , k 2 E )), 2 ^ with pe the momentum in 
Euclidian space. The divergence of g{p 2 E ) appearing at pe = ^qcd, is removed by introducing 
an infrared cutoff parameter piF 2 ^' 2 ^ as 

f{vl) = ^((pI + pIv^cd) ' (24) 



with 

1 1 1 AT _ 9 /V, 

(2-5) 



1 11JV C -2JV> 



8tt 2 3 

Here, iV c and iVj are the numbers of colors and flavors, respectively. 

We next introduce the following bilocal auxiliary field for non-perturbative analysis: 

Xapip, q)=J K a ^ 5 (p, k)^ (k + |) j s (k - |) . (2-6) 

This makes the action bilinear in the quark fields if) and if). Then, integrating out if) and if), 
we obtain the classical action Z — J T>xexpi r[x] as 

r lx] = \ I ^[x(P,~(l)K~ 1 {p,k)x(k,p)] 

^ J pkq 

-*Tr Ln(# 4 (g)(27r) 4 - X (p, -?))■ (2-7) 



The SD equation is obtained as the extremum condition for this classical action within 
the stationary-phase approximation for x- 13 ^ The existence of a non-trivial solution for x 
indicates the dynamical breaking of the chiral symmetry. 

Due to the translational invariance of the vacuum, this solution is expressed as (x<xp) — 
^'a / 9(p)5 4 (9)(2vr) 4 , with the mass function E, and the effective potential V[E] = — r[ X ]/ J d 4 x 
becomes 



V[E] = -\ ( tr[E{p)K-\p,k)E{k))+i [ tr \n(^-E(p)). 

^ J pk Jp 

Then the SD equation, SV/SE = 0, in the improved ladder approximation is written 



(2- 
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Whereas E a p{jp) has the general form E{jp 2 )8 a p + E v {p 2 )^ a g in the vacuum, the S v {p 2 ) part 
is known to vanish in the Landau gauge (a = 0). 27 - 1 After carrying out the Wick rotation 
and the angle integration, the SD equation becomes 

where C 2 (N C ) = T a T a = (N 2 - 1)/2N C . 

2.2. Order parameters of the chiral symmetry 

The quark condensate with the four momentum cutoff A is defined as 

This bare value at the scale A is converted into the value at the lower energy scale fi (e.g., 
1 GeV) via the renormalization group equation 



i 

AD 



W>m = W>a(^J , (2-12) 

where B = (12C 2 (N C ))- 1 (11N C - 2N f )/3. We will present the value of (^)m at /i = 1 GeV 
and omit the subscript \x hereafter. 

Next, the pion decay constant f n is estimated in terms of the mass function U(p 2 ) by 
utilizing the Pagels-Stokar formula: 28 ^ 

N c f™ p 2 E E(p 2 E ) / , 2 . p\dE{p\) 



n = M, «ra7f«-^ ^ (2 ' 13) 

We fix the value of Aqqd here so as to reproduce the empirical value of /„. with this formula. 

2.3. At finite temperature and chemical potential 

We use the imaginary time formalism to extend the QCD-like theory to the case with 
finite temperature and chemical potential, making the following replacement: *) 

/ f{Po,P) ► T J2 [ **f(iu n + p,p), (2-14) 

P n=—oo 

where ui n = (2n + l)nT (n G Z) is the Matsubara frequency for the fermion. 

*) We ignore the possible dependence of the running coupling constant g 2 on the scales, T and/or ^i. 29 ) 
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The mass function E a p(u n ,p) at finite T and //, is invariant under spatial 0(3) rotations 
and decomposes into E(co n , \p\)S a p + E s (u n , |p|)o> n (7o)a/3 + E v (u n , \p\)p l {^i) a p. It is in fact 
possible, 20 ^ although still cumbersome, to solve the SD equation numerically in this general 
form. For our purpose of demonstrating the usefulness of the effective potential in the QCD- 
like theory, we assume here E s = E v = for simplicity. We, further, use a covariant-like 
ansatz 26 ) for the mass function, taking 



E(uj n ,p) — > E(p 2 ), 



(2-15) 



where the frequency and the momentum appear in the combination p 2 = cj 2 + \p\ 2 in E. 
Thus, the SD equation simplifies to 



E(f 



3C 2 (N C 



8tt 5 



- g 2 (P 2 ,k 2 ) 



dk l 



ra= — oo m 



2 — w 2 



xln 



p 2 + k 2 + 2\ (p 2 - w 2 n )(k 2 - w^) - 2w n w r 



x 



p 2 + k 2 — 2\J (p 2 — w 2 l )(k 2 — — 2w n w r , 
k 2 + 2ifiu m - fi 2 + E{k 2 ) 2 ' 



(2-16) 



At finite chemical potential, solutions of this equation are generally complex-valued. We set 
lu h = on the right-hand side of Eq. (j2-16|) . which makes the SD equation real. At high 
temperature (T ^> \p\ and T ^> |fe|), the Matsubara frequencies with large n give only small 
corrections to physical quantities. 2&> 

Substituting the solution E into the effective potential (|2-8|) with the replacement (|2-14j) . 
the extremum value Vl x reads 



n=0 " U n 
2 ,.2 



(p 2 -^ 2 + E 2 (p 2 ))E 2 (p 2 ) 
(p 2 - fx 2 + E 2 (p)) 2 + V^ r 2 



\n[(p 2 -^ 2 + E 2 (p 2 )) 2 + 4^ 



(2-17) 



The stability of the symmetry-broken phase is usually examined by comparing this extremum 
value with that of the trivial solution. It is preferable to have a functional form of the 
effective potential in order to study the features of the phase transition. In the next section, 
we present a method for determining the shape of the effective potential for the composite 
field. 
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§3. Shape of the effective potential 



3.1. The effective potential away from the extremum 

Here we explain a method that we use to construct the functional form of the effective 
potential numerically. A standard way to assess the potential form is to apply an external 
source which is coupled to the field linearly. It is recognized, however, that using this 
approach we cannot study non-convex potentials, which is an important feature near the 
phase transition point. 

We apply an external source field J(p, k) which is coupled to the square of the self-energy 

as 

V[E, J] = V[E] + ~ / tx[E(p)J(p, k)E(k)}, (3-1) 

Z Jpk 

where V[E] has been defined in (|2-8jl . By imposing the extremum condition for V[E, J] with 
respect to E(p) (that is, 5V[E, J]/8E{p) = 0), we derive the SD equation with the source 
field J as 

We denote the solution of ()3-2j) by Ej(p) to indicate the implicit dependence on the source 
J. The effective potential V for the configuration Ej(jp) is written 

V[Sj) = V[Ej, J] - \ [ U[Ej{p)J{p, k)Ej{k)} 

= ~lJ p tT Sj ^ \-E % j(p) +i J p tTln ^- S j(p)) 

-\ [ tr[Ej(p)J(p,k)Ej(k)}. (3-3) 

z Jpk 

We study the global behavior of V[Zj] near the critical points by varying Ej through a 
particular type of variation of a function J(p, k). 

Here, among infinite possibilities, we consider one natural choice, J(p, k) = —cK~ l (p, k), 
with a parameter c, and study the shape of the effective potential along this particular 
variation. Then, the SD equation becomes 

£e(p) = yj- I K(p, k) 1 (34) 
l + cj k fk- E c (k) 

We note that the variation of c is, in effect, equivalent to varying the strength of the strong 
coupling constant, since (K^ 1 — J) -1 = K/(l + c) oc g 2 /(l + c). In the limit c — ► oo we 
should have a trivial solution, while a symmetry-breaking solution, E c 0, exists in the 
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case c — > — 1, which corresponds to the case of infinite effective coupling. In this way, 
we should be able to scan the potential shape along a family of configurations, from a trivial 
one to a symmetry-breaking one. 

By substituting the solution E c {p) of ()3-4j) into V[S C ) in (jH Hjl . we obtain the effective 
potential function of c: 



7(c) 



21 + c. 



/ tr Z c (p) I + i [ trlnd/ -E c (p)). (3-5) 



As the scalar condensate (ipip) is computed with U c , the effective potential can be expressed 
numerically as a function of (V>V)- Near a first-order transition point, the SD equation (J3-4)) 
has two non-trivial solutions for a certain range of c, each of which corresponds to a different 
value of the scalar condensate. Using both the solutions, we obtain the effective potential 
for the full range of the scalar condensate. It is straightforward to extend this potential to 
a system at finite temperature and chemical potential using the replacement expressed in 

(Ell- 
in the numerical evaluation of the effective potential, we split V[E] as 

V[Z) = {V[i:)-V[0))+V[0). (3-6) 

Here, the difference between the free energies V[Zj] — V[0] is finite and can be evaluated 
numerically. V [0] = i f tr In $ is the potential of a free massless quark gas and is divergent 
due to the vacuum fluctuations. It is thus necessary to adopt an appropriate regularization 
to remove this divergence. In our approximation, we replace the last term V[0] in ()3-6|) with 
the pressure Pf ree of a free massless quark gas: 



■V[0] - P free = N c N f T 4 



7tt 2 + 1/^\2 + 1 fl i\4 



180 6 \TJ 12tt 2 \T 



(3-7) 



In studying the nature of the phase transition, especially near a tri-critical point, it would 
be instructive to extend the effective potential to a form with two independent order param- 
eters, the quark number density and the scalar density. 31 )***-* To this end, we first replace 
the variable \i with the quark number density p through the usual Legendre transformation, 

F((#),p,T) = V{(W),H,T)+pti , (3-8) 

*' This situation is similar to that in the strong coupling QED investigated previously. 30 ' 
**) For this particular choice for J, one could also utilize the Hellmann-Feynman theorem. We acknowl- 
edge referee's comment on this point. 

***> In the case of the QCD critical point with a finite current mass, this extension of the potential becomes 
more useful. 
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ln(p 2 /AL) 



Fig. 1. The mass functions £(p) as functions of p. The left panel displays £(p) with p = MeV 
at T = 0, 50, 80, 100, 110 and 115 MeV, respectively. The right panel displays £(p) with 
T = 80 MeV at p = 0, 280, 290 and 300 MeV, respectively. 



with 

p = -^((^),//,T) . (3-9) 
a p 

Then, by adding a coupling energy with an external potential v to Eq. ()3-8|) . we define a 
new effective potential as 

V((W), p; v, T) = F((W),P, T) - pv % (3-10) 

which may be interpreted as a Landau potential. 32 > In order to construct a non-convex 
function of p we use the unstable solution of the SD equation. Note that p and v are 
independent variables here. When the condition dV /dp = p — v = is imposed, this new 
potential V coincides with the original potential: V((i(j'^), p(p); p, T) = V ((tpip) ; p,T). 

3.2. Numerical results 

We studied the case with two massless quark flavors (i.e., Nf = 2) and N c = 3, where 
N c is the number of colors. The model parameter Aqcb is fixed so as to reproduce the 
pion decay constant /„. = 93 MeV in ()2-13|) . and we obtain Aq CD ~ 776 MeV using the IR 
cutoff at ln(pj F /ylQ CD ) = 0.1. Then, the quark condensate in the vacuum is found to be 
| (?/>?/>) I 1 / 3 ~ 236 MeV in this model, and it is insensitive to the value of the UV cutoff. (We 
use A 2 /A 2 QCD = 2.0 x 10 4 .) 

At finite temperature and density, we numerically calculated the quark mass function 
E (p) and the quark condensate (ipip). In Fig. E]we plot the mass function S(p). The left 
panel displays S{p) with p = MeV at T = 0, 50, 80, 100, 110 and 115 MeV, respectively, 
while the right panel displays E(p) with fixed T = 80 MeV at p = 0, 280, 290 and 300 MeV, 
respectively. At the finite chemical potentials p = 280 and 290 MeV with T = 80 MeV, 
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Fig. 2. The chemical potential (left) and temperature (right) dependences of ^(0) and Ktpip)^^ 3 , 
respectively. The left figures plot these quantities at T = 50,80, 100 and 110 MeV as functions 
of jJL. The right figures plot these two quantities at ji = 0, 200 and 300 MeV as functions of T. 



there exist two non-trivial solutions for E(p), which we denote I and II. The free energy of 
the solution I is lower than that of II. 

In the left panels of Fig. |2J we plot the \i dependences of £"(0) and {ipip), respectively, 
at T = 50, 80, 100 and 110 MeV. The right panels present the T dependences of the same 
quantities at /i = 0, 200 and 300 MeV. In the region of low temperature and large chemical 
potential, we find two non-trivial extremum solutions, £ ^ 0, in addition to the trivial one, 
which indicates a first-order transition. It should be noted here that the chiral condensate 
(ipip) and the value of E(0) increase as /i increases in the case of fixed T = 50 MeV, as 
is seen in Fig. |21 This behavior is known to be an artifact of the approximation made by 
omitting E s and E v . If E s and E v are properly taken into account, (tptp) and E(0) are found 
to decrease monotonically, and the chiral symmetry is restored, as /x increases. 20 ) 

We present the phase diagram of the model in FigEl where the tri-critical point appears 
at (T t ,/i t ) ~ (97,203) MeV. 

The situation may be more clearly understood if we can obtain the explicit functional 
form of the effective potential. We plot V[E C ] — V[0] as a function of Aq CD in Fig. El 

at /j, = 0, 100 and 300 MeV for several values of T. In Fig. we plot the results at T = 80 
and 100 MeV with several values of \x. As seen from the behavior of the effective potential 
in Fig. HJ a second-order transition occurs in the range T = 100 - 120 MeV in the cases with 
\i = and \i = 100 MeV. At /i = 300 MeV, however, a first-order transition occurs between 
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600 



H [MeV] 



Fig. 3. The phase diagram for the chiral symmetry. The transition point is located at T c ~ 114 
MeV for y, = MeV, \i c ~ 472 MeV at T = MeV, and so on. The tri-critical point P is 
located at (T t ,Ht) ^ (97,203) MeV. 




Fig. 4. The temperature and chemical potential dependence of the effective potential, V[E C , J] 
V[0], with fixed \i as a function of 



T = 75 MeV and 80 MeV. Similarly, in Fig. El a first-order transition is seen to occur when 
the chemical potential is changed between fi = 280 MeV and 290 MeV at low temperature 
(T = 80 MeV), while the transition is second order at fixed T = 100 MeV. In Fig. the 
quark number density p is plotted in units of the normal nuclear matter density in terms of 
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Fig. 5. The temperature and chemical potential dependence of the effective potential, V[27 C ) J] 
V[0], with fixed T function of \(ipip}\. 




(i [MeV] 

Fig. 6. The chemical potential dependence of the quark number density at T = 80 and 100 MeV. 
Here, p = 3x 0.17 fm~ 3 . 



quark numbers, po — 3 x 0.17 fm 3 , as a function of p with fixed T = 80 MeV and 100 MeV. 
The transition points are p = 282 MeV at T = 80 MeV and p = 187 MeV at T = 100 MeV. 
In the former case, the quark number density has a gap. In the latter, it is continuous but 
its derivative (the susceptibility) has a gap. 

In Fig. we plot the contour map of the Landau potential V((ipip), p; p, T) given in 
(|3T0J) as a function of the quark condensate and the quark number density at several values 
of p with fixed T. Note that a tri-critical point is located at (Tf,p t ) ~ (97,203) MeV. The 
vertical axis represents (ipip) / Aq CD , and the horizontal axis represents the quark number 
density in units of po- We plot the contour map at T = 100 MeV, where a second-order 
transition occurs, in the left panels of Fig|7| At p = 0, the two minima are positioned 
symmetrically in the p-(ipip) plane at p = with finite values of the quark condensate, (^pip). 
As the chemical potential increases, the two minima approach each other for finite p, and 
then fuse continuously to form a single minimum at (ipip) = and finite p. 

The contour map for T = 80 MeV is plotted in the right panels of FigO At p = 0, the 
two minima are positioned symmetrically in the p-(ipip) plane at p = and finite values of 
the quark condensate, (ipip)-, just as in the T = 100 MeV case. By contrast, however, one 
new local minimum appears at (ipip) = as p increases, and three local minima exist in a 
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Fig. 7. The contour map for V({ijjijj),p ; /i, T) in the p-fyi/j) plane with various chemical potentials 
for T = 100 MeV (left) and T = 80 MeV (right). 



certain range of /x. At the critical point, these three minima are energetically degenerate, 
and we observe a first-order chiral phase transition. 

§4. Summary 

We have analyzed the chiral phase transition of a QCD-like theory at finite temperature 
and density using the effective potential. We have devised a method to derive an effective 
potential in order to elucidate the global behavior of the phase transition. There we in- 
troduced a bilocal external source field and solved the Schwinger- Dyson equation with this 
source field. This solution of the SD equation gives the extremum of the effective potential 
with the source field. Then, subtracting the interaction energy with the source we obtained 
the original effective potential for this configuration. In this way we showed the shape of 
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the effective potential of the QCD-like theory at finite temperature, T, and quark chem- 
ical potential, /i. From the temperature dependence of the effective potential, we have a 
second-order phase transition along the n = line. Contrastingly, we have a first-order 
phase transition along the T = line. Also, we have introduced the Landau potential and 
plotted its contour map with respect to the quark condensate and the quark number density. 
We conclude that our method to construct the effective potential of the QCD-like theory is 
very useful for understanding the features of the chiral phase transition of the model. 
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Appendix A 

Calculation of the Effective Potential in the NJL Model 



For pedagogical reasons, here we apply our method with a source field coupled quadrat- 
ically to the auxiliary field, to the Nambu-Jona-Lasinio (NJL) model, whose Lagrangian 
density takes the form 

£ = ^ + A ((^) 2 + (fer^) 2 ) • (A-l) 

We introduce the auxiliary field into the scalar part of the 4-fermi interaction and ignore the 
pseudo-scalar part for the moment, as we know that the potential is chirally symmetric. To 
leading order in the 1/N expansion, we obtain the effective potential at finite T and fi as 



y 2 f d :i k 

V(T, ft, X ) =^ - 2N c NfT / 7^T-3 ln [K - + k 2 + X 2 

n=—oo ' 



y 2 n N f r x r A k 2 

X c f ' dx d\k\ (l-n+(x) -n.(aQ), (A-2) 



2A 2vr 2 J J JW^~x 



where x = ~ Xipip/N is an auxiliary field, and n±(x) = [exp(\/ k 2 + x =F /i) + l] -1 . 

The functional form of the potential is easily obtained through the numerical integration 
of Eq. (|A-2[) . To demonstrate our method, we consider the case of a constant source coupled 
to the field y 2 ; 33 ) 

V(x,J) = V(x) + \jx 2 - (A-3) 
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Fig. 8. Two effective potentials V(x) — V(x = 0) (solid curve) and V(x, J) — V(x 
(round points) composed in the case T = 10 MeV and = 320 MeV. 

Rewriting the source J as J = cA -1 , we obtain 



0, J = 0) 



V( X , J) = (1 + c)f- - -%-± / dx d|fc|-^=(l - n + {x) - n_(*)), (A-4) 
2A 2n 2 ./o Jo Jk 2 + x 



and the corresponding SD equation, dV /dx = 0: 



Xc 



X 1 + c 



1 iV c iV /Xc 



7T" 



d|fe| 



fc 2 



(l- n+ (xl)-n4x 2 c ))- 



(A-5) 



The value of the effective potential at x — Xc is obtained from V ex (xc) using the Legendre 
transformation: 

1 



v(xc) = v cx ( X c) - ~Jxl ■ 

In the stationary-phase approximation for x, we should have the trivial result 

v( x ) = v( x ). 



(AS) 



(A-7) 



One of the advantages of the source (jA-3|) is that it allows us to probe the non-convex 
effective potential near a second-order transition. Near a first-order transition point, the 
SD equation (jA-5|) has two non-trivial solutions for a certain range of c. We need both the 
solutions in order to re-construct the full shape of the potential via Eq. (jA-6|) . 

We present the result in the massless two-flavor case (A^ c = 3 and Nf = 2) here. We fixed 
the three-momentum cutoff as A ~ 653 MeV and the coupling constant as XA 2 ~ 4.29 with 
the constituent quark mass 313 MeV and /V = 93 MeV. In FigJHl we compare the effective 
potential V(x) (circles), constructed using the applied source and the Legendre transforma- 
tion, with the original V(x) (solid curve). It is seen that they are essentially coincident, 
as should be the case. Hence, we conclude that our method successfully reproduces the 
potential in the non-convex region. 33 ) 
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